A new contribution to the nuclear modification factor of 
non-photonic electrons in Au+Au collisions at s/snn = 200 GeV 



G. Martmez-Garci'a^*, S. Gadrat^''' and P. Crochet^* 
^SUBATECH (IN2P3/CNRS - Ecole des Mines - Universite de Nantes) Nantes, France 
Laboratoire de Physique Corpusculaire (IN2P3/CNRS - Universite Blaise Pascal) 

Clermont-Ferrand, France 
(February 2, 2008) 

Abstract 

We investigate the effect of the so-called anomalous baryon/meson enhance- 
ment to the nuclear modification factor of non-photonic electrons in Au+Au 
collisions at ^snn = 200 GeV. It is demonstrated that an enhancement of 
the charm baryon/meson ratio, as it is observed for non-strange and strange 
hadrons, can be responsible for part of the amplitude of the nuclear modifica- 
tion factor of non-photonic electrons. About half of the measured suppression 
of non-photonic electrons in the 2 — 4 GeV/c pt range can be explained by 
a charm baryon/meson enhancement of 5. This contribution to the non- 
photonic electron nuclear modification factor has nothing to do with heavy 
quark energy loss. 
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One of the most robust experimental evidence for the creation of a new state of matter in 
heavy-ion coUisions at the Relativistic Heavy Ion CoUider (RHIC) is the large suppression of 
light hadrons at high transverse momentum (pt) [!]• This phenomenon is well reproduced by 
models which take into account the radiative energy loss of high pt light quarks and gluons 
propagating through a dense medium of colored quarks and gluons [2] . Further insights into 
the underlying mechanism can be obtained from the study of heavy hadrons. Indeed, heavy 
hadrons originate from the fragmentation of (heavy) quarks in contrast to light hadrons 
which are predominantly produced by gluons. Quarks are supposed to lose less energy 
than gluons in the medium due to a smaller color charge coupling. In addition, radiative 
energy loss was predicted to be substantially smaller for heavy quarks as compared to light 
quarks because of the so called "dead-cone" effect which limits the medium induced radiative 
energy loss at forward angles [3] . Surprisingly, recent data from the PHENIX and the STAR 
collaborations in Au-|-Au coUisions at ^/sNN — 200 GeV show that the quenching of non- 
photonic electrons lies beyond theoretical expectations [4-6] and is as large as that of light 
mesons. This is difficult to understand under the common assumption that non-photonic 
electrons originate predominantly from the decay of D and B mesons, the latter originating 
themselves from heavy (c and b) quarks hadronization. Reconciling these data with model 
predictions is a real challenge which triggers a lot of theoretical activities nowadays. Only 
models which assume a very large medium opacity [7] or an additional coUisional energy 
loss [8,9] or a very early fragmentation of heavy quarks in the medium [10] can describe 
successfully the data (for a recent review, see [11]). It is worth noticing that, in contrast to 
light hadrons, the heavy flavor quenching is, so far, not measured experimentally through 
identified particles, but in an inclusive way via the nuclear modification factor (Raa) of 
non-photonic electrons^ . 

In this paper, we point out the possibility that part of the explanation for the strong 



The nuclear modification factor of non-photonic electrons is obtained from the pt distributions 
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suppression of non-photonic electrons might be due to a another source of electrons namely 
charmed baryons. Indeed, whereas light mesons are largely suppressed in heavy ion collisions 
at RHIC, the suppression of non-strange and strange baryons is observed to be much less 
in the intermediate pt range (2 < pt < 4 GeV/c) [12]. This is commonly referred to as 
the anomalous baryon/meson enhancement. It could challenge our understanding of jet 
quenching at a first glance but brings however another piece of evidence for the formation of 
a deconfined medium. This anomalous baryon/meson enhancement is well understood in the 
framework of the recombination model which assumes that, at moderate pt, hadronization 
occurs via the coalescence of "free" quarks (and anti-quarks) [13-18]. In the following, we 
show that an anomalous baryon/meson enhancement for charm hadrons leads naturally to 
a non-photonic electron Raa smaller than one. This is essentially due to a smaller semi- 
leptonic decay branching ratio and to a softer decay lepton spectrum of charm baryons 
as compared to charm mesons. As a consequence, a fraction of the deviation from the 
experimentally measured Raa of non-photonic electrons should not be attributed to energy 
loss. After the description of the physical process invoked here, we present simple PYTHIA- 
based simulations to illustrate the magnitude of the effect in Au-|-Au collisions at -s/snn — 
200 GeV and we compare our results to the PHENIX data. 

The main assumption we put forward is that, in a deconfined medium, charm baryon pro- 
duction is enhanced relative to charm meson production, as compared to the vacuum. This 
assumption is qualitatively justified in the framework of the recombination model. Although 
this model does not provide predictions on charm hadron production yet, it successfully de- 
scribes the baryon/meson enhancement measured in Au-|-Au collisions at y^Sjvjv — 200 GeV. 



of non-photonic electrons in AA collisions (dA^^^/dpt) and in pp collisions {dNpp/dpt) as: 

dN%Jdpt 
< N^J > dNe^/dp, 

where < > is the average number of nucleon-nucleon collisions corresponding to a given 

centrality class. 



A relatively good agreement is obtained not only for the light hadron ratio p/'K'^ i but also 
for heavier hadron ratios such as A/X^ and Vt/cf) [21]. However, extrapolating these results 
to charm hadrons is not straightforward because the mass of the charm quark is much larger 
than that of hght and strange quarks. The consequences are threefold as far as the recom- 
bination mechanism is concerned. First, whereas the px of a light baryon (meson) amounts 
to 3 (2) times the initial pt of its valence quarks, the pt of a (single) charm baryon or a 
charm meson is likely to be very close to that of the charm quark. Secondly, considering 
a light quark and a heavy quark with the same velocity (which is the essential require- 
ment for the coalescence process to take place [22]), the heavy quark momentum is much 
larger than that of light partons. As a consequence, one can expect the enhancement of 
charm baryon/meson to appear at higher pf The recombination model indeed predicts, 
for non-strange and strange hadrons, that the heavier the hadron, the larger the pi of the 
baryon/meson enhancement [21]. Finally, the fragmentation time of heavy quarks is small 
as compared to light quarks. According to [10], the formation time of a 10 GeV/c pion, D 
meson and B meson is 20, 1.5 and 0.4 fm/c, respectively and it is as small as ~ 3 fm/c 
for a Ac with = 20 — 30 GeV/c. Due to these considerations, it is obvious that the 
baryon/meson enhancement for non-charm hadrons and charm hadrons can be significantly 
different. We stress that our aim is neither to investigate these differences nor to quantify 
the magnitude of the charm baryon/meson enhancement. This would require a sophisti- 
cated theoretical investigation which goes clearly beyond the scope of this paper and, to 
our knowledge, no information on this topic is available in the literature yet. Therefore, in 
the following, we only assume that, in view of experimental results on the baryon/meson 
enhancement for non-strange and strange hadrons, a similar enhancement is a priori con- 
ceivable for charm hadrons. Remarkably, such an enhancement has strong implications on 
the nuclear modification factor of non-photonic electrons. It leads to a decrease of the yield 
of non-photonic electrons in AA collisions because, as shown in Tab. I, the inclusive semi- 
leptonic decay branching ratio of charm baryons is smaller than that of charm mesons (and 
the hadronization process does not change the total charm quark number). Therefore, the 
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nuclear modification factor of non-photonic electrons should decrease as well. This can be 
easily illustrated in the following way. Assuming that a AA collision is a perfect sum of pp 
collisions (i.e. Raa = 1) and that the relative yields of D mesons are the same in pp and in 
AA coUisions, apt integrated Raa can be calculated for different Ac/D enhancement factors 

C = N^^^^JNd = {Na^ + N^J/{Nd+ + Nd- + Nno + Nj^o + Nj,+ + Nj,-) according to 

1 + C{Nj,^,j^jNn) 1 + iVA,^e/iVz,^e ^ ' 

where 

AT /AT ^ i^Ac,Ac/^D)BRA,,A, 

{Nd±/Nd)BRd± + {N^o^do/Nd)BR^o^d^ + {N^±/Nd)BR^± ' ^ ' 
N is the charm hadron yield and BR is the hadron semi-leptonic decay branching ra- 
tio. According to Tab. I, N^^^/^JNd = 7.3%, Nd±/Nd = 21%, Nj^o^do/Nd = 67%, and 
Nj^±/Nd = 12% such that NA,-.e/ND-.e = 3.63%. Enhancement factors C of 5 and 12 
therefore lead to a Raa of 0.90 and 0.79 respectively. It is shown below that the effect is 
further amplified in the intermediate pt range as a result of a softer decay electron spectrum 
of Ac as compared to that of D mesons. 

In order to investigate this effect in more details, we have performed simple simulations 
with the PYTHIA-6.152 [23] event generator. The PYTHIA input parameters were first 
tuned according to [20] and the PHENIX acceptance cut (|?7| < 0.35) was apphed in order to 
correctly reproduce the pt distribution of non-photonic electrons measured in pp collisions at 
^/s = 200 GeV [24]. As it can be seen from Fig. 1, the agreement between the simulation and 
the data is rather good except at very low pt and in the high pt region where the simulation 
under-predicts the data. The disagreement in the high pt region is due to the fact that 
we consider, in the simulation, only electrons from charm decay and neglect electrons from 
bottom decay which are expected to dominate the spectrum for pt ~ 4 GeV/c [24]. On 
the other hand, it can be observed from Tab. I that the Ac/ D ratio amounts to 7.3% (in 
47r) and, as mentionned above, to 3.63% after convolution of the species yields with their 
corresponding semi-leptonic decay branching ratio. Figure 2 shows that, in the 2 — 4 GeV/c 



Pt region of interest discussed hereafter, this ratio is even smaller (~ 1.5%) because the decay 
electron spectrum of Ac is softer than that of D mesons. In this respect, an enhancement 
of the Ac/ D ratio in central heavy ion collisions induces two different mechanisms of non- 
photonic electron yield suppression: i) integrated yields are suppressed due to the lower 
semi-electronic branching ratio of the Ac and ii) the softer electron pt distribution from Ac 
decay leads to an additional suppression at intermediate pt- 

We then re-evaluate the non-photonic electron decay spectrum after introducing a Ac/ D 
enhancement chosen to 5 and 12. This enhancement is assumed to be fiat in pt and it is 
applied such that the pt-differential charm cross-section is conserved. The latter is an arbi- 
trary choice that could be justified since most of the charm hadron transverse momentum 
is given by the charm quark whatever, baryon or meson, this hadron is. We finally compute 
the Raa ratio from the previous non-photonic electron p^ spectra neglecting all effects like 
shadowing or quenching but assuming that the only medium induced effect is the Ac/ D 
enhancement. The nuclear modification factor is therefore simply obtained from the ratio 
of the non-photonic electron pi distribution with a Ac/ D enhancement to the original distri- 
bution. By doing so, one implicitly takes into account the pt dependence of the AJ D ratio 
which is not fiat (Fig. 2). The results are shown in Fig. 3 together with the PHENIX data. 
The simulated Raa ratio is shown only for 2 < pt < 4 GeV/c since shadowing and electrons 
from bottom decay, which play an important role at lower and higher pt respectively, are 
not considered in the simulation. One can see from Fig. 3 that a A^/ D enhancement of 5 
can already explain ~ 50% of the suppression of non-photonic electrons in the 2 — 4 GeV/c 
Pt range. For illustration, we show that, in the extreme scenario of an enhancement of 12, 
the data can be satisfactorily described without invoking any heavy-quark energy loss. The 
shape of the simulated Raa is obviously determined by the pt distribution of Ac/ D enhance- 
ment that we have assumed to be flat. Note that the values of Raa at pt ~ 3 GeV/c can be 
roughly obtained as well from Fig. 2 and Eq. 1 and 2, assuming that the decay electrons at 
Pt ~ 3 GeV/c originate from charm hadrons with Pt ^ 4 GeV/c. 

In summary, we have shown that an enhancement of the Ac/D ratio in heavy ion coUi- 
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sions has dramatic consequences on the nuclear modification factor of non-photonic electrons. 
Such an enhancement, which is observed for non-strange and strange hadrons, would signif- 
icantly lower the Raa of non-photonic electrons (or muons) at intermediate Pf This results 
from i) the semi-leptonic decay branching ratio of charm baryons which is smaller than that 
of charm mesons and ii) the decay lepton spectrum from charm baryons which is softer than 
that of charm mesons. We conclude that it is therefore premature to interpret the non- 
photonic electron Raa data before a possible enhancement of the Ac/ D ratio is measured 
experimentally and/or investigated theoretically. Heavy quark energy loss can be studied in 
a much cleaner way from the nuclear modification factor of exclusively reconstructed charm 
hadrons. Such measurements should be possible with the RHIC-II experiments [25] and 
with the ALICE experiment at the LHC [26] . We finally note that the Ac/ D enhancement 
can possibly influence the elliptic flow of non-photonic electrons as well. 

NOTE 

After the publication of the present work in arXiv we have realized that a similar work 
has already been published by P. Sorensen and X. Dong in PRC 74 (2006) 024902. However, 
the assumptions made in the two approaches are different and lead to different conclusions. 
In PRC 74 (2006) 024902, assuming a AJD ratio of 1.7, similar to the measured A/K° 
ratio, about 35% of the of non-photonic electron Raa data is reproduced. In our approach, 
the Raa data are reproduced by using the Ac and D spectra from PYTHIA and assuming 
a AJD ratio of ~ 1 (i.e 12 x 0.08). 
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TABLES 

TABLE I. Inclusive decay branching ratio (BR) of charm hadrons into e + anything [19] and 
yield (A) of charm hadrons (in An) in pp collisions at -^s = 200 GeV from the present PYTHIA 
simulation using input parameters as described in [20] . The total cross-section for charm production 
is normalized to the experimental value obtained in [20] . Aa^ and A^^ include primarily produced 

Ac and Ag as well as those from and Eg decay. 



Hadron 


D+ D- 


DO DO 


Dt 




Ac Ac 


BR (%) 


17.2 ± 1.9 


6.71 ± 0.29 






4.5 ± 1.7 


A (xlO-^) 


3.00 3.07 


9.31 9.85 


1.82 


1.60 


1.23 0.85 
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FIGURES 




(GeV/c) 

FIG. 1. Invariant differential cross-section of non-photonic electrons (dots) measured in pp 
collisions at ^/s = 200 GeV [24]. The solid curve shows the result of the PYTHIA simulation as 
described in the text. The simulated spectrum is normalized from the integration of the measured 
spectrum in the range 1.4 < pt < 4 GeV/c. 
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FIG. 2. Charm baxyons over charm mesons (dots) and decay electrons from charm baryons over 
decay electrons from charm mesons (squares) versus pt , from the PYTHIA simulation described in 
the text. 
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FIG. 3. Nuclear modification factor of non-photonic electrons (dots) measured in central 
(0-10%) Au+Au collisions at ^/snn = 200 GeV [4]. The solid and dotted curves correspond 
to the results of the simulation described in the text for a Ac/D enhancement factor of 5 and 12 
respectively. 
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